Ge-substitutional defects and the ^/S x y/3 ^-> 3x3 transition in a— Sn/Ge(lll) 
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The structure and energetics of Ge substitutional defects on the a-Sn/Ge(lll) surface are analyzed 
using Density Functional Theory (DFT) molecular dynamics (MD) simulations. An isolated Ge 
defect induces a very local distortion of the 3x3 reconstruction, confined to a significant downwards 
displacement (-0.31 A) at the defect site and a modest upward displacement (0.05 A) of the three Sn 
nearest neighbours with partially occupied dangling bonds. Dynamical fluctuations between the two 
degenerate ground states yield the six-fold symmetry observed around a defect in the experiments 
at room temperature. Defect-defect interactions are controlled by the energetics of the deformation 
of the 3x3 structure: They are negligible for defects on the honeycomb lattice and quite large for a 
third defect on the hexagonal lattice, explaining the low temperature defect ordering. 
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The driving force behind the T-induced reversible 
y/3 X y/3 <-> 3 X 3 phase transition observed for 1/3- 
ML coverage of Sn or Pb on Ge(lll) has been under in- 
tense debate since its discovery in Pb/Ge(lll) |l). After 
the early suggestions of surface Fermi wavevector nest- 
ing or electron correlations, leading to the formation of 
a Charge Density Wave ||l|;p| at low T, recent theoretical 
work indicates that the softening of a 3x3 surface phonon 
, associated with the electronic energy gain due to the 
Sn dangling bond (DB) rehybridization, plays a major 
role on this transition. The dynamical fluctuations model 
g associated with this soft mode provides a consistent 
account of the main experimental facts, in particular the 
similarities and differences between the isovalent Sn/Ge 
and Sn/Si(lll) systems |]. 

DFT calculations for the a-Sn/Gc(lll) surface pre- 
dict that the ground state is a 3x3 reconstruction in 
which 1/3 of the Sn adatoms are placed ~ 0.3 A higher 
than the other Sn adatoms |^-^. In other words, out of 
the three 3x3 sublattices that form the \/3 x v^ lattice 
| fiO[ , "up" Sn adatoms (Sriu) occupy one of these sub- 
lattices (say a), forming a 3x3 hexagonal lattice, while 
sublattices 6, c are occupied by 'down' Sn adatoms (Srid) 
(forming a honeycomb pattern). These calculations also 
show that this reconstruction is correlated with a change 
in the surface electronic structure H : the Sn^-DB is now 
fully occupied while the other two S'nd-DBs share the re- 
maining electron per 3x3 unit cell. STM experiments are 
naturally explained in terms of these different DB occu- 
pancies. At LT, below Tc, filled-state STM images show 
the 3x3 phase : bright spots occupy the 3x3 sublattice 
associated with the S'n,(-DBs, while the other two sub- 
lattices, associated with the S'^^-DBs, are occupied by 
darker spots 0. At RT, however, dynamical fluctuations 
of the Sn adatom heights, fast with respect to the STM 
scan speed, induce a \/3 x a/3 periodicity in the STM 
images. 

While these models assume an ideal 1/3-ML Sn sur- 



face coverage, point defects, with an average concentra- 
tion of 3±1 % Q, are always present. 90 % of these 
defects are Ge substitutionals (called Ge defects in the 
rest of the paper), while vacancies account for the other 
10 %. STM, a real-space experimental technique, has 
been the method of choice to get non-averaged informa- 
tion on the structure and distribution of these defects on 
Sn/Ge(lll) [pl^p^ across a wide temperature range, cov- 
ering the transition temperature Tc- Three are the main 
experimental results: (i) At RT Ge defects are imaged 
as dark spots in the v^ x \/3 lattice surrounded by six 
features (brighter than Sn adatoms away from defects) 
that correspond to the six nearest-neighbours (n.n.) Sn 
atoms, (ii) At LT most of the Ge defects occupy positions 
corresponding to the two honeycomb S'rid-sublattices 6, c, 
with equal probability, and are therefore surrounded by 
three bright and three darker Sn n.n. (iii) At RT Ge de- 
fects are randomly distributed over the three sublattices 
a, 6, c. 

These results have been considered as supporting evi- 
dence for an interpretation of the \/3 x v3 ^^3x3 phase 
transition in terms of the interaction of the Ge defects 
|&-nO[. One of the main ingredients in this argument 
is the six-fold symmetry assumed for the perturbation 
induced by an isolated defect. This six- fold symmetry 
(honeycomb, with one atom dark and two atoms bright in 
the three possible sublattices), suggested by the RT STM 
images, is, however, in contradiction with the hexagonal 
pattern (two atoms dark and one atom bright, with three- 
fold symmetry) observed in the filled-state images below 
Tc- In order to explain the LT experimental results, a 
defect-defect interaction, with a temperature dependent 
decay length 1{T) (e.g. 1{300K) ~ 11 A, 1{100K) ~ 100 
A) has to be introduced. This long-range interaction 
would be responsible for the ordering of defects at LT on 
two of the three possible sublattices. Finally, the super- 
position of the perturbations induced by defects in these 
two sublattices would give rise to the 3x3 phase. 



In this letter, we show that ab initio DFT-LDA 
MD simulations of Ge substitutional defects on the 
Sn/Ge(lll) surface support a completely different point 
of view: Both the structure and the ordering of the Ge 
defects at LT are induced by the ground state 3x3 re- 
construction. In particular, we have studied the atomic 
geometry, the electronic structure and the energetics of 
one and two Ge defects in a large 9x9 unit cell. In the 
case of two defects, we have considered three cases, cor- 
responding to the defects on first, second and third n.n. 
positions. Our analysis shows that Ge defects at LT al- 
ways create a three-fold distortion around them, the sim- 
ulated STM filled-states image showing three bright and 
three darker spots on the n.n. Sn atoms. The distor- 
tion is very local, affecting only to the Sn atoms which 
are nearest neighbours of the defects. We have also per- 
formed first-principles MD simulations of a 3x3 unit cell 
with one Ge defect. Our MD simulations indicate that 
the six-fold pattern observed around Ge defects in the 
RT STM images should be understood in terms of the 
RT dynamical fluctuations of the six n.n. Sn adatoms. 
These dynamical effects contribute, together with struc- 
tural and electronic effects, to the enhanced brightness 
of the Sn atoms around the defect. 

In our approach we calculate different defects using a 
Sn/Ge(lll)-9x9 unit cell with a slab of 4 Ge layers and 
a layer of H atoms saturating the dangling bonds of the 
deeper Ge layer (432 atoms in the unit cell). Both the 
H layer and the last Ge layer are kept fixed during the 
calculations to mimic the bulk continuation. In our cal- 
culations we have used an efRcient MD technique, the 



Fireball96 code |12|, that allows us to analyze the large 
9x9 unit cell we have considered. We have checked that 
this code yields accurate surface geometries by compar- 
ing it with a plane wave code (CASTEP) ||ll[ for the case 
of the Sn/Ge(lll)-3x 3 reconstruction 1^ and for the case 
of one Ge defect on a 3 x 3 surface unit cell (see Table I) . 
These results confirm the quality of the calculations pro- 
vided by the Fireball96 code and justify the wide use we 
are going to make of it in the rest of the paper. STM im- 
ages are calculated from the ab-initio electronic structure 
of the surface using the Keldysh-Green function formal- 
ism and considering a W tip (for details on the method 
and the application to the Sn/Ge(lll)-3x3 system see 
ref. §). 

Figure la shows the relaxed geometry obtained in our 
calculations for the a-Sn/Ge(lll) surface with one Ge 
defect per 9x9 unit cell (3.7 % concentration). The sur- 
face presents the 3x3 reconstruction with the Ge adatom 
on one of the two 3 x 3-sublattices (say b) that form the 
honeycomb pattern associated with Sn^ adatoms. Our 
simulations yield for the Gc substitutional a downwards 
displacement of 0.30 A with respect to the original Srid 
position. The perturbation induced by the Ge defect 
is well-localized, affecting only to the three n.n. Sud 
adatoms (c sublattice). These adatoms are displaced up- 



wards ^ 0.05 A due to the transfer of charge from the 
Ge defect: the Ge-DB is empty and each n.n. S'rid-DB 
gains ^ 1/6 electrons. The occupancy of the n.n. 5'n„- 
DBs (a sublattice) does not change (they are already fully 
occupied), and therefore their position is practically the 
same as in the 3x3 reconstruction free of defects. Fig- 
ure Qp shows the calculated filled-state STM image for 
this geometry. Probably the most remarkable feature is 
the three-fold symmetry of the image around the defect 
(only three of the surrounding Sn atoms are imaged as 
protrusions), in contrast with the six- fold symmetry as- 
sumed for the perturbation induced by the defect from 
the experimental RT STM images. 




FIG. 1. (a) Relaxed structure for the isolated defect. Only 
the vertical displacements (in A) with respect to the ideal 
3x3 structure which are larger than 0.01 A are indicated. 
Notice that the distortion is confined to the defect site (black 
circle) and the three neighbouring Sn atoms with partially 
occupied dangling bonds (gray circles). Sn atoms with fully 
occupied DBs are indicated by white circles, (b) Filled-state 
STM image for the isolated defect. The bright protrusions 
correspond to the fully occupied Sn DBs, and the gray ones 
to the Sn atoms with partially occupied DBs. The image 
reflects the three- fold symmetry of the ground state structure 
of the defect. 

Notice that the solution shown in fig. [^ is two-fold de- 
generate: an equivalent solution is obtained applying a 
symmetry operation that exchanges the two sublattices 
(a, c) occupied by the up and down Sn. This degeneracy 
will be important when discussing the effect of tempera- 



ture on the surface geometry and its corresponding STM 
image. 

Figures 2a,b,c display the relaxed adatom positions we 
have obtained for the a-Sn/Ge(lll)-9x9 surface with 
two defects on first, second and third n.n. positions. In 
all the cases the 3x3 reconstruction adapts its particu- 
lar geometry (out of a three-fold degenerate case) to the 
defect positions, in such a way that the two defects are 
placed on the two S'rxd-sublattices forming a honeycomb 
pattern. When the Ge defects are second n.n. (fig. 2b) 
they are on the same sublattice (say 6), and the ground 
state is two-fold degenerate: an exchange of the a and 
c sublattices gives a completely equivalent solution. In 
fig. 2a and 2c the Ge defects occupy different sublattices 
and the solution is non-degenerate. As in the case of an 
isolated defect, the perturbation induced by the defects 
is localized essentially on the Sn-adatoms n.n. to the 
defects. The total energy for the three cases shown in 
fig. g is, within the accuracy of these calculations (~ 20 
meV/9x9 unit cell), the same, confirming the localized 
nature of the distortion induced by the Ge defects. It is 
remarkable that the sum of the upward displacements of 
all the neighbouring Sn-adatoms is, in the three cases, 
very similar to the sum of the displacements induced by 
two independent defects (~ 0.30 A, see fig. pi). This 
result suggests a roughly linear relationship between the 
upwards displacement of a Sn adatom and the transfer 
of charge to its DB. 

Coming back to the experimental STM data for the Ge 
defects on the Sn/Ge(lll) surface, we first notice that 
our solution for a single defect (fig, lb) has a three-fold 
and not a six-fold symmetry, with the Sn atoms around 
the defect presenting two different types of displacements 
perpendicular to the surface. This apparent contradiction 
with the experimental data at RT can be resolved, how- 
ever, considering the dynamics of the atoms around the 
defect. To this end, we have considered a 3x3 unit cell 
with a single substitutional and performed MD simula- 
tions at two different temperatures, T= 180 K and T= 
270 K. The initial conditions correspond to the positions 
of the relaxed structure of the defect (see Table I) and 
random velocities following a Maxwell-Boltzmann distri- 
bution. Then we follow the time evolution of the system 
integrating the equations of motion with the forces ob- 
tained from the DFT-LDA calculations. Figure || shows 
the time evolution of the z-coordinate of the three upper 
atoms (two Sn atoms and the Ge susbstitutional). These 
results show that Ge remains around its equilibrium site, 
while the two Sn atoms tend to exchange heights, very 
much in similarity with what we found for the \/3 x v3 
structure M. Therefore, the RT STM images of the six 
Sn-adatoms n.n. to a Ge defect must represent a 50 
% average of Sriu and Snd (displaced upwards r^ 0.05 
A w.r.t. normal S'/irf-adatoms, and containing 1/6 elec- 
trons more), while the RT STM images of Sn-adatoms 
away from defects represent a 1/3 Sriu + 2/3 Sud av- 



erage. This dynamical (plus charge transfer) effect ex- 
plains in a simple way the different pattern observed at 
RT and LT on the STM images around a Ge defect (six- 
fold vs. three-fold symmetry) and also explains why Sn- 
adatoms n.n. to a Ge defect appear on the RT STM 
images brighter than Sn-adatoms in regions free of de- 
fects. 




FIG. 2. Relaxed structures (as in fig. hi) for two Ge de- 
fects on first (a), second (b) and third (c) nearest neighbour 
positions. 

We finally address the following question: why do Ge 
defects at LT occupy preferentially two 3x3-sublattices, 
with equal probability, and not the three of them or only 
one?. Our results show that, on the a-Sn/Ge(lll) sur- 
face, Ge substitutional defects transfer their DB charge to 



neighbouring Sn-adatonis DBs and do not participate in 
the up/down dynamical fluctuations of the Sn-adatoms 
at RT. Our results also show that at LT, when these 
fluctuations are supressed, the stable configuration of the 
surface is a 3x3 reconstruction in which Ge defects oc- 
cupy Sud positions; the defect-defect interaction is neg- 
ligible (as long as they stay on Sud sublattices) and Ge 
defects occupy with equal probability any of the two Sud 
sublattices. 
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FIG. 3. Time-evolution of the 2;-coordinate of the surface 
atoms (two Sn atoms and the Ge susbtitutional -dotted line-) 
on the 3x3 unit cell in MD simulations at 180 and 270 K. 

A third Ge defect initially located in the other sub- 
lattice (hexagonal) should be forced, when the temper- 
ature is lowered, either to present an upward displace- 
ment (contrary to what it does when it is isolated) or 
to deform dramatically its environment in order to move 
downwards. Any of these solutions would imply a sig- 
nificant energy cost (160 meV for the defect in an upper 
position according to our calculations), that the system 
prefers not to pay, locking in all the substitutionals in 
a honeycomb lattice. We conclude that the ordering of 
Ge defects at LT is induced by the ground state 3x3 
reconstruction and not the other way around: Ge substi- 
tutional defects are not the driving force of the VS x VS 



^^3x3 transition in the Q:-Sn/Ge(lll) surface. 

In conclusion, our DFT-LDA total-energy and MD cal- 
culations show that the ground state geometry around 
a Ge substitutional defect corresponds to a very local 
distortion of the 3x3 reconstruction of the defect-free 
surface, confined to the defect site and the Sn nearest 
neighbours. The lattice deformation can be character- 
ized by a downwards displacement of the Ge defect and 
upward displacement of the neighbouring Sn atoms with 
partially occupied DBs, with respect to the ideal 3x3 
structure. This structural effect is correlated with the 
charge transfer from the defect DB to the DBs of these 
Sn atoms. This ground state configuration provides, at 
LT, an STM image with a three-fold symmetry, while 
at high T, the dynamical fluctuations between the two 
degenerate states (with different Sn atoms in the upper 
positions) would yield the six-fold symmetry observed 
in the experiments. These dynamical effects, together 
with the upward displacement and the increased DB oc- 
cupancy, are responsible of the enhanced brightness of 
the Sn atoms around the defect. Finally, the ordering of 
Ge defects at LT is naturally explained as an effect of the 
3x3 reconstruction. 
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TABLE I. Atomic vertical displacements of the adatoms 
(in A), measured with respect to the 3x3 reconstruction, for 
a Ge defect in a 3x3 unit cell. For the Ge defect the dis- 
placement is referred to the position of one of the "down" Sn 
atoms of the 3x3 reconstruction. 



Sn "up" 



Sn "down" 



Ge defect 



PW 
F96 



+0.04 
+0.02 



+0.08 
+0.10 



-0.31 
-0.32 



